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An approach to the calculation of intermodulation noise in FM systems 
due to imperfect transmission media is presented in this paper. The tech- 
nique is essentially that originating with Carson and Fry. In this paper we 
extend a formulation due to Rice to include an arbitrary continuous pre-em- 
phasis characteristic as well as an arbitrary gain and phase shape transmis- 
sion medium which are representable by low-order polynomial series in radian 
frequency. Series approximations are carried out far enough to ensure ac- 
curate results for transmission characteristics normally encountered in 
broadband microwave radio systems. In many cases, only the second- and 
third-order noise is significant in broadband microwave radio systems. Hence, 
the analysis carried out in this paper considers only the second- and third- 
order distortion terms. A digital computer program concerning the inter- 
modulation noise has been written. This analysis and the digital computer 
program are of aid in the design of microwave 1 adio systems. With a slight 
modification, the calculation of noise due toAM-w-PM conversion caused by 
transmission deviation can also be accomplished. An optimum design of 
the pre-emphasis network may be achieved by using the computer programs 
through an iterative approach. 

I. INTRODUCTION 

In the course of designing FM systems, ir.termodulation noise is an 
important factor which deserves special attention. Many people have 
made contributions to this subject. 1 - 10 In this paper, we extend their 
results to include an arbitrary continuous pre-emphasis characteristic 
as well as an arbitrary gain and phase shape transmission medium which 
are representable by low-order polynomial series in radian frequency. 
Series approximations are carried out far enough to ensure accurate 
results for transmission characteristics normally encountered in broad- 
band microwave radio systems. The multichannel baseband signal of an 
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FM system is represented by a band of gaussian random noise with flat 
power-density spectrum. The noise due to imperfect transmission me- 
dium can be calculated at any frequency in the baseband. In many 
cases, only the second- and third-order noise is significant in broadband 
microwave radio systems. Hence, the analysis carried out in this paper 
considers only the second- and third-order distortion terms. Extension 
to a higher-order distortion becomes unmanageable. A digital computer 
program concerning the intermodulation noise has been written. A 
typical problem can be solved at a very low cost. This analysis and the 
digital computer program are of aid in the design of microwave radio 
systems. With a slight modification, the calculation of noise due to 
AM-to-PM conversion caused by transmission deviation can also be 
accomplished. An optimum design of the pre-emphasis network may be 
achieved by using the computer programs through an iterative approach. 
Several examples are given for illustration. 

II. DESCRIPTION OF SYSTEM 

A portion of an FM system can be represented by the block diagram 
shown in Fig. 1. An FM baseband signal, <p/ (t) = d<pi(t)/dt, is fed into 
a pre-emphasis network. The output of the pre-emphasis network is 
the pre-emphasized signal <p it). This signal passes through an FM 
modulator to yield the FM wave, cos [w c t + v(t)] where co c is the angular 
carrier frequency. When this FM wave goes through an imperfect 
transmission medium, the output is distorted and becomes V(t) cos 
[o) c t + <po(t)]. Let the transfer function of the transmission medium be 



F(co) = exp [-«(/) -10(f)], 
where <a = 2tt/. Its impulse response is 

g (t) = I Y(u) exp (iut) df. 



(1) 



(2) 
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BASEBAND SIGNAL 
<P'(t) 



FM WAVE 
COS [(U C t + 9>(t)] 
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V(t)cos[wct+<Po(t)] 



TRANSMISSION 
MEDIUM 




Fig. 1 — Block diagram of a portion of an FM system. 
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Using the complex notation, the input and output of the transmission 
medium can be written, respectively, as 

vt(t) = exp [iw c t + *p(0] (3) 

and 

vo(t) = 7(0 exp [iio c t + fo(t)], (4) 

where 7(0 is taken to be positive and <po(t) is determined to within 
2nr, n is an integer. 

The output and input of the transmission medium are related by the 
convolution integral 

vo(t) = I Vi(t - x) g{x) dx. (5) 

J— oo 

From (3), (4), and (5) 

7(0 exp [i<po(t)] = I exp [vp (t — x) — iio c x] g(x) dx. (6) 

J — 00 

Let 

7(0 = exp[a(0]. (7) 

Then 

a(0 = Re In / exp [i<p(t — x) — io) c x] g(x) dx (8) 

J— oo 

<po(t) = Im In / exp [i<p(t — x) — iu c x] g(x) dx. (9) 

J— 00 

The AM distortion term expressed in dB is 20 X 0.4343 a(i) and the 
PM distortion term expressed in radians (or degrees) ispo(£) — <p(t)- 

Assume that the transmission medium passes only frequencies in 
the neighborhood of the carrier frequency, ±/ c ± b, with b/f c « 1. 
Thus, to a high degree of approximation, we have 

exp (— iu3 c x) g(x) ~ I Y(u c + w) exp (itax) df. (10) 



Let 



k(x) = ^-. — r exp (— iu c x) g(x). (11) 



Substituting (1) and (11) into (8) and (9) we obtain 
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a(t) = -oc(f c ) + Re In / exp [up(t - x)] k(x) dx (12) 

J— 00 

<p {t) = -/3(/ c ) + Im In / exp [up(t - x)] k(x) dx. (13) 

J— 00 

The quantity k(x) may be regarded as the normalized envelope func- 
tion of the impulse response g(x). We shall discuss now the logarithm 
of the integral in (12) and (13) in detail. 

III. DERIVATION OF DISTORTION TERMS 

Let 

exp [vp{t - x)] = M(t,x). (14) 

A delay t d is often introduced in order to improve the degree of approxi- 
mation of various series with the first few terms. Then M (t,x) can be 
expanded about x = td as 



M(t,x) - M(t,t d ) + {X 1 , td) \^M(t,x)\ 

1 ! \_dx Jx=t 



[£ m «'4 



+ tUt* m «-* ) .L + 



(15) 



One choice of td is 



t d = Re £ xk(x) dx = ± [*$"L • ( 16 ) 

Using (14) and (15), the integral in (12) and (13) can thus be expressed 



as 



f exp M* - x)] k(x)dx = £ ™f f~|- M(t,x)~\ , (17) 
where m n is the nth. moment of k (x ) defined as 

m n = I (x — td) n k(x)dx, wo = 1 (18) 

J— 00 

or equivalently, 

fc^T \^-Tn r(«c + «) exp (f«ft,)l . (19) 

r(w c ; \_d{zo}) n J u =o 



w„ = 



The series (17) is equivalent to the Carson-Fry series with delay t d . 
It may not converge in certain cases. However, when the characteristic 
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of a transmission medium can be truly represented by a polynomial, 
from (19), the higher moments become zero and the series reduces to a 
polynomial. The logarithm of (17) can be written as 



In / exp[i<p(t — x)] k(x)dx = i<p(t — td) 

J — en 



|_ «=i n\ J 



(20) 



where 



F n (t - t d ) = (-1)" exp [-Mt - t d )\ [J^te*)] • 
Using Taylor's series expression, we have* 
In (l + S «n(t)x n /nl\ = j- { ai + fj («■ - «i 2 > 

3 

+ o] («3 — 3aia 2 + 2ai ) 

4 

+ — (a 4 — 4aia 3 — 3a 2 2 + 12«i 2 a 2 — 6ai 4 ) + • • • . 
4! 

With x = — 1 and a n (t) — m n F n (t — td), substituting the above series 
into (20), after considerable algebra, one obtains an asymptotic series 
which has been written by S. O. Rice in an unpublished work as 

f° r- /. m,/ \. T m\ i . mi a mz in 

In J exp [up(t - x)\ k(x)dx = i\<P ~ ±1* + 2\ * ~ $[ * 

■ mi mi i i \ '3 

+ XT" * —•"•+* W 

— % (m 4 — 2wim3 —mi 

+ 2m 1 2 m 2 W'+---l (21) 

+ [-iX#> /2 + |(m 3 

— mimi)<p'<p" — |(m 4 

— mim 3 )<p(p — f (m 4 

— mi)<p"~ + ^* X 4 ^' 4 + • • • ], 

* Notice that this is not the approximation In (1 + j/) ^ y which can be a poor 
approximation and yet has been quite widely used in FM work. 
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where <p stands for <p(t — t d ) and \„'s are the semi-invariants which 
are related to the moments by 

X 2 = m 2 — mi 

X 3 = m 3 — 3mim 2 + 2m] 3 

\i = mi — 4rai7tt 3 — 3ra 2 2 + 12wi 2 m 2 — 6rai . 

Taking the real and imaginary parts of (21 ) and substituting the result 
into (12) and (13) gives, respectively, 

a(t) = - a{f e ) + m U tp — -^f +"37^ — "JT * + * * ' 

X3i /3 , hi 12 II X 2r n , Z 2r / /' /nn\ 

— — ^ + -j<p*> — ~2 <P -r-^vv \**) 

hr I HI Ur 112 1 A4 r n . 

-j<P<P -j<P + ^<P + ••• 

*»(0 ■■ ~ P(f*) + v ~ «iif +"2]"*' ""31"^ "^"iT*' -••• 

. X3r '3 Zlr / 2 // X 2 f / 2 i 2 ,' / // hi I HI /oO\ 

~j<P +24^ + •••, 

where the subscripts r and i denote the real and imaginary parts of the 
corresponding coefficients and 

li = ra 4 — 2miW 3 — m% + 2rai 2 ra 2 , 

h = ma — WI1WI2 , 

k = nii — mim 3 , 

k = rtii — ra 2 2 . 

From (22) and (23), amplitude and phase distortions are divided into 
linear, second- and third-order terms and are shown in Table I. The 
terms —a(f c ) and —p(f c ) in (22) and (23) do not appear in Table I. 
This is due to the fact that they are constants and introduce only con- 
stant amounts of amplitude and phase distortions. 

IV. PRE-EMPHASIS CHARACTERISTIC 

The multichannel baseband signal of an FM system is represented 
by a band of random noise. Assuming that the bottom baseband fre- 
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Table I — Amplitude and Phase Distortions Due to 
Imperfect Transmission Medium 



Order of Distortion 


Amplitude Distortion a(<) 


Linear 

Second-order 

Third-order 


mu<P - — <? + — <? -—<? +■■■ 

X2r , 2 hr , „ hr , „, hr „, 

- T f' 3 + 7 *V' + • • ■ 

6 4 




Phase Distortion ¥*)(/) — <p(t) 


Linear 

Second-order 

Third-order 


--*■ + -«»* - e** " g* + ••• 

\ir , hr , „ 

— <p * — -r v *<t> + ■•• 
6 4 



quency is much smaller than the top baseband frequency, the power- 
density spectrum of the baseband FM signal is expressed as 

where /& is the top baseband frequency. 

A pre-emphasis network is used in an FM system in order to optimize 
the noise across the baseband. Let Z(a>) be the transfer function of the 
pre-emphasis network, we write 

| Z(u) | ' = a + a 2 / 2 + a 4 f + a 6 / 6 , | / | ^ fb , 

where the a's are real constants either given a priori or determined by 
the least squares fitting from an actual curve. The power-density spec- 
trum of the pre-emphasized baseband FM signal is 

S v .(w) = |£(o>)l 2 <V(«) 

= P (oo + <hf + a 4 f* + a 6 f), \f\ ^ f b . 

In the case when the pre-emphasis coefficients a , a 2 , 04 , and o 6 
are determined by the least squares fitting from an actual curve, the 
following weighting function of normalized scale has been found useful 
for better approximations near the bottom channels 



(24) 
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W(f/f b ) = io-'""™"^- 11 , 

where Z is the difference of relative power (in dB) of top and bottom 
channels of the given pre-emphasis characteristic. 

The rms frequency deviation, <r, due to noise loading can be expressed 
as 

(27r«r) 2 = ave [<p ,2 (t)] = [ S f >(a) df. 

•J — oo 

Using (24), the relation between P and a can be expressed as 

■Po = j= / Tax / 7T\ / ,6\-i , (rad/sec) 2 /Hz, 

where the units of <r and /& are Hz. 

V. TRANSMISSION MEDIUM 

Within the band of interest, f c ± b, let the gain and phase of the 
transmission medium be, respectively, 

exp [-«(/ + f e )] = 1 + giu + 92U + g*o? + g& 

» (25) 

+ z2 Uk cos (p k u + 0*) 

iV 

-£(/ + / c ) = &2co 2 + b 3 J + 6 4 co 4 + Zft sin (#a> + (r fc ), (26) 

where the gf's and &'s represent coarse shape transmission deviation; 
the u's and v's represent fine shape transmission deviation. It should be 
emphasized that the fine shape transmission deviation is restricted to 
be a slowly varying ripple characteristic, hence this analysis does not 
apply to the noise due to single echo in general. The fine shape representa- 
tion is useful when we study the effect on the noise of a given system 
due to the shift of the carrier. Since a constant delay does not introduce 
intermodulation noise, the linear term in co is not included in (26). The 
transmission medium coefficients g, b, u, v, p, q, 6, and a in (25) and (26) 
are either given directly or obtained approximately by a curve-fitting 
computer program. 

Substituting (26) into (16), we have 

w 
td = — JL M* cos *k • 
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From (19), we can evaluate various moments, hence, the coefficients 
associated with the distortion terms in Table I in terms of transmission 
medium coefficients, as given in Appendix A. In the following sections 
we shall derive expressions for intermodulation noise calculation due to 
second- and third-order distortion terms. 



VI. NOISE POWER DUE TO SECOND-ORDER DISTORTION TERM 

Since 

(I '- o ./ " 

-tp = liptp , 

at 



cU- 



d n n i in | o in 

— tp = Zipip + Zip , 



the second-order PM distortion in Table I can be written approximately 



as 



where 



»(*) = (-*x« + \hi j t - Mi ;j-j <p' 2 + AW 2 , (27) 



la — 4/sf 3t5t . 



The second-order PM distortion term of (27) can be represented by the 
block diagram shown in Fig. 2, where 

Hi(w) = (x&ZsiW 2 — 1*2*) + *(&«»), 

//o(o)) = A^4* • 

The power-density spectrum of ^ 2 (0 is 1 

S M (w) = Fi(-o>)Hi(w)S/i(«) + ffi(-«)ifr(«)jS,'V'«(«) 



4- fft(-«)Hi («)iS f »■/■(«) + H 2 (-u)H2(w)S p "i(a), 



(28) 



H.(w) 



fz(t) 



H,(w) 



Fig. 2 — Block diagram representation of the second-order PM distortion term. 
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where S^tita) and S v »2(u>) are the power-density spectra of <p' 2 and 
<p" 2 , respectively, and S V '2 V »2(u) [or S 1 ,»i r >*(a)] is the cross-power- 
density spectrum of <p'~ and ip" 2 (or ip" 2 and <p' 2 ). These quantities can be 
derived as 12 

SM") = W(r) +72/(0)], 

flf^V'«(«) = 3 ! PBr'r» f (r) + B^ (0)fl,» (0)], 

S,»V.(w) = ff[2fl,. v 2 (r) + /2,'(0)fl,''(0)], 

£„".(«) = JF[2^- 2 (r) + Z2„» 2 (0)], 

where R v -{t) and R v >>(t) are the autocorrelation functions of <o' and 
<p", respectively, and R V ' V "(t) [or ijVv'(r)] is the cross-correlation func- 
tion of <p' and <p" (or p* and <o')> and 5? stands for "the Fourier transform 
of". 

The Fourier transform of a constant function is a delta function at 
zero frequency. In the situation of evaluating noise power in the base- 
band, this quantity is not of interest. Also, it can be shown that 

S v '* v "a(o)) = 8 r "a v 'i(<a). 

Thus, (28) can be simplified as 

fl w (») - 2 | Hi{») | 2 *[V(r)] + 2 | ffi(«) | '^[R v -\r)\ 

+ 2[ffi(-o>)# 2 (co) + Ft(-«)ffi(«)MB^» f (r)]. 

The intermodulation noise to signal power ratio due to the second order 
distortion term expressed in dB is, therefore, 

N 2 /S(dB) = 10 log [JS Vi (u)/S 9 . («)]. (30) 

VII. NOISE POWER DUE TO THIRD-ORDER DISTORTION TERM 

Since 

- v =3<p<p , 
the third-order PM distortion in Table I can be written approximately as 

V*(t) = [&3r — Tshr 37 J ^ ■ 

The above equation can be represented by the block diagram shown in 
Fig. 3, where 
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H 3 (w) 



Fig. 3 — Block diagram representation of the third-order PM distortion term. 

H 3 (u) = (|X 3 r) + t(-tWl^). 

The power-density spectrum of p 8 (0 is 

S V3 (<») = I #3(w) | 2 S v '*(u), 

where *Sys(w) is the power-density spectrum of <p' 3 which can be derived 
as 12 

£„'s(w) = 6JF[#/(r)] + 9«/(0)^(co). 

In the above equation, the term 9R V > (0)>Sv (co) is merely a scaled power- 
density spectrum of the input baseband FM signal, hence, it does not 
contribute to the intermodulation noise and can be neglected in the 
computation. 

The intermodulation noise to signal power ratio due to the third-order 
distortion term expressed in dB is, therefore, 

N 3 /S(dB) = 10 log [coX ( a )/S, («)], (31 ) 

where 

S,» -6|JI,(«)|'ff[V(T)]. (32) 

In (29) and (32), the Fourier transforms of R v > (r), R v >- (t), R v > v >> (t), 
and R v - 3 (t) may be obtained by taking the convolutions in the frequency 
domain. However, this requires numerical integration. For given pre- 
emphasis characteristics and P (or a), these Fourier transforms can be 
expressed in algebraic forms as shown in Appendix B. Hence, no nu- 
merical integration is necessary. A digital computer program has been 
written to calculate the second- and third-order intermodulation noise 
due to second- and third-order distortion terms in dB. A typical problem 
can be solved at a very low cost. 

VIII. EXAMPLES 

Several examples are considered in this paper. Calculated results are 
compared with measured data when they are available. Expressions for 
noise calculation are derived for simple cases. For more complicated 
situations, the noise calculation is best carried out by using a digital 
computer. 
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8.1 Example 1 

In this example, we wish to demonstrate how the intermodulation 
noise across the baseband can be optimized using an appropriate pre- 
emphasis network. To simplify the calculation, we assume that the 
transmission characteristic consists of linear delay distortion (6 2 ) only. 
From Appendix A, we obtain 

X 2t - = -2b 2 ,hr = -86 2 2 

and all the other coefficients are equal to zero. Hence, 

ffi(«) = b 2 , HM = 0, H 3 (a) = if&2 2 «. 

In practical cases, the third-order distortion term [/f 3 (w)] is negligible 
(say, 50 dB less) compared with the second-order distortion term. From 
(30) we write 

* ( dB ) = 10 log Willi . 

o o V '(w) 

For simplicity, we let 

Mco) = P (l + a,/ 2 ), |/| £/*. 
From Appendix B, we obtain 

S[iV(r)] = Po*f„v(f), 



where 



Since 



v(f) = —hAW - f A<$? + (fill + 2)A&? 
- (1 + A 2 ) 2 fl + (tW + U2 + 2) 
A 2 = ctifb 2 
0-//A. 



p ^ (2«r) 



4+4-°)' 



consequently, 



^(dB) = 10 log 



(2r)W»T 7 



©W,) 



+ t) 



^ki + Arfi") 
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Specifically, we let 

bo = 7.962 X 10 -17 (linear delay of 1 nanosec/MHz) 

f b = 1 MHz 

a = 1 MHz. 

After several computer runs, the optimal choice of a 2 is 7. The inter- 
modulation noise with no pre-emphasis and with the optimal pre-empha- 
sis are plotted in Fig. 4. The noise has been reduced to more than 5 dB 
and is evenly distributed across the baseband by using the optimal pre- 
emphasis network. 
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Fig. 4 — Intermodulation noise due to linear delay with and without pre- 
emphasis. 

8.2 Example 2 

In this example, we Use a typical radio system pre-emphasis charac- 
teristic shown in Fig. 5. The top baseband frequency is fb = 5.772 MHz. 
Since the given pre-emphasis characteristic is expressed as relative 
power in dB versus baseband frequency, it is first converted to ratio 
versus normalized baseband frequency, ft = f/fb ■ 

The weighting function is 

W(Q) = ior* w-K,w) . 

A least squares approximation program is used to obtain the approxi- 
mating polynomial 

no + d-4- + «4.f 4 + a 6 f 6 , 
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where 

a = 0.99894166 

a 2 = 11.944252// 6 2 

a, = -5.5771705// b 4 

a 6 = 1.4396088//6 6 . 

Consider a transmission characteristic consisting of a linear, a parabolic, 
and a slowly varying sinusoidal delay as shown in Fig. 6. The expres- 
sions for the noise due to second- and third-order distortion are too com- 
plicated to write down. However, by using a digital computer, the 
results are plotted in Fig. 7 for a = 0.771 MHz. Clearly, a better 
pre-emphasis network should be used to optimize the noise across the 
base-band for this particular transmission characteristic. 

8.3 Example 3 

As a final example, we consider a single pole IF filter with 

1 



F((0 + co c ) = 



w 



where w is the 3-dB half-bandwidth of the filter. Using a least squares 
approximation with appropriate weighting function, the magnitude and 
phase of F(oj + u c ) are expressed by 
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Fin 5.— Pre-emphasis characteristic of a typical radio system. 
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Fig. 6 — An arbitrary delay characteristic of a transmission medium. 
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Fig. 7 — Intermodulation noise due to the delay distortion of Fig. 6. 



1552 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1966 

exp [-«(/ + /.)] ^ 1 + g*£ + g¥ ** = 1 + &(//«>)» + G5,(f/«)« 

-/3(/ + /„) £* & 1U + 6 3 co 3 = Bi(Vw) + # 3 (t» 3 , 

where 

G 2 = j/ 2 (27riy) 2 = -0.4209, G 4 = g^w) 1 = 0.08027 

j?x = b!(2irM>) = -0.9529, B 3 = b 3 (2irw) 3 = 0.1294. 

The actual and approximated transmission characteristic are plotted in 
Fig. 8. Since a constant delay does not introduce intermodulation noise, 
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Fig. 8 — Gain and phase characteristic of a single pole filter. 
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61 is not included for calculation. Using the expressions in Appendix A, 
the coefficients associated with the second-order distortion term are 



Hence, 



X 2 ,- = 0, hi = 0, l 3i = 0, la = 0. 



ffi(w) = // 2 (a.) = 0. 



The noise contribution due to second-order distortion is, therefore, zero. 
The coefficients associated with the third-order distortion term are 

X 3 r = 663 , kr = 2404 - 40 2 2 . 

Assuming no pre-emphasis, that is, a 2 = 0,4 = a 6 = 0, from Appendix 
B, we have 



where 



Using the relation 



$[K' 3 (r)] = Po 3 / 6 W(3 - ft 2 ), | ft | ^ 1 
ft = ///* • 
(2tt<t) 2 



P a = 



2/ 6 a 



the inteiTnodulation noise to signal power ratio due to the third-order 
distortion term can be derived from (31) as 

*(dB> = 101ogl© 4 g)V (3 - rf) [b; + (£)'(«, - f )]. 

For/b = 1 MHz and w = 1.25 MHz, N 3 /S (dB) is calculated at/ = 0.084, 
0.36, and 1 MHz as a function of (o// 6 ). The dotted lines in Fig. 9 repre- 
sent the calculated value (S/N s ) while the solid curves represent the 
measured data taken by W. F. Bodtmann. 13 The discrepancy between 
the measured and calculated values can be attributed to several reasons: 
(i) The power-density spectrum of the multichannel baseband signal 
used in the experiment was not perfectly rectangular, (it) During the 
measurement, a non-ideal limiter was used which caused some AM-to- 
PM conversion, (Hi) The actual transmission characteristic was approxi- 
mated by few parameters in a limited region, and (iv) the formulas (30) 
and (31) derived in this paper involved approximations. Nevertheless, 
the measured and calculated results are close enough to show the utility 
of the analysis even for cases beyond the application for which it was 
originally intended. 
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Fig. 9 — Measured and calculated intermodulatiou noise due to the single 
pole filter. 



IX. DISTORTION DUE TO AM-TO-PM CONVERSION 

Let K be the AM-PM conversion factor of a device expressed in de- 
grees/dB, then the phase distortion due to AM-PM conversion is 

Phase Distortion due to AM-PM Conversion = 20 X 0.4343K a(t) 

= 8.686it~ a(t) degrees, 

where a(t) is given by (22). 
Similarly, 

Frequency Distortion due to AM-PM Conversion = -^ — K — ■ a(t) 

2iir dt 

= 1.382 K 4 o(0 Hz. 

dt 

In Table I, the second- and third-order terms of amplitude distortion 
cause intermodulation noise while the linear term causes video roll-off 
or enhancement. Using the same approach discussed previously in this 
paper, noise calculation due to AM-to-PM conversion which is caused 
by transmission deviation can also be accomplished. 



NOISE IN FM SYSTEMS 1555 



X. ACKNOWLEDGMENTS 



The analysis carried out in this paper is based on an unpublished work 
of S. O. Rice. His permission to present his result here is greatly appreci- 
ated. The author is grateful to J. T. Bangert for his encouragement and 
direction and also wishes to thank J. J. Lang, F. J. Witt, L. H. Enloe, 
C. L. Ruthroff, and W. T. Bamett for many constructive suggestions. 
The author is indebted to Miss J. D. Witkowski for the digital computer 
programming and to many other colleagues who have helped to carry 
out this analysis. 

APPENDIX A 

Coefficients Associated with the Distortion Terms 
The moments defined in (19) can be expressed as 
Wo = 1 



•G-) 



nil =t r 



™ 2 = (-^j + *(-*") 






nu 
where 

N 

A = l + E«* cos fa 



A' = 0i — X) u k pic sin d k 

fc=l 

N 

A" = 2#2 — 2 u k -p k cos e k 
fc=l 

N 

A'" = 6^3 + S u k p k 3 sin d k 



A-=l 



A" ' = 2404 + J2 u k p k cos 6 k 
fc=i 
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B" = 2b 2 — 2 v k qk sin <r k 

N 

B'" = 66 3 — 2 v k qk cos <r k 

N 

B"" = 24& 4 + Z mn sin «r t . 



The coefficients associated with the second-order distortion are 
X21 = — B 

A'" , i'l' 

D ,„, , 3A'B'" , QA"B" ZA n B" 
La = B H j -f- 



A ' A A 2 

12A"£" <B iV 



fo = B"" + ^7^- - 12 



A A 2 

The coefficients associated with the third-order distortion are 

X 3r = B'" 

A"" a" 2 a' a'" a' 2 a" 



A A 2 A 2 A 



APPENDIX B 



#Wier Transforms of R v - 2 (t), fl, V 'V), R v » 2 (t) and R v > 3 (t) 

The Fourier transforms of R r - 2 (r), flpvV), #„" 2 (t) and iV( T ) can 
be derived in a straightforward manner. However, considerable amount 
of algebra has been involved during the derivation. Without writing out 
the details here we merely present the final results. 

S[R v > 2 (t)] = Po%{D 1 (f - 2f b ) + 2D 1 (f) + D 1 (f + 2/0 

- JhU - 2/») + 2D 2 (f) - D 2 (f + 2/ 6 ) 
+ 2D,C/-2/») -2Z) 3 (/+2/ 6 )| 

5[R..*(t)] = Po 3 f b 2 {SE 1 (f - f b ) + Ei(f - 3/ 6 ) + jEitf + 3/ 6 ) 
+ SBtf + /a) + 3# 4 (/ - fb) -82E»(f - 3/ 6 ) 

- 3£ 4 (/ + 3/ 6 ) + 3# 4 (/ + /») + 3tf ,(/ - ft) 
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+ 3E 3 (f - 3/ b ) 

- ZE 3 (f + 3/ 6 ) - 3E 3 (f + ft) + 3E*(f - f b ) 

- E % (j - 3/ 6 ) + £»(/ + 3/ b ) - 3E& + /„)} 
JF[i2„ V ' 2 (r)] = 47r 2 P 2 / b 3 { Jj(/- 2/ b ) + 2/i(/) + Ji(/+ 2/ b ) - J 2 (/- 2/ b ) 

+ 2J 2 (f) ~ U1+ 2/ fc ) - 2J 3 r.f - 2/ 6 ) + 2J 3 (/+ 2/ b )} 
S[*V' 2 (t)] = 16tt 4 PoWKiif - 2/ b ) + 2tf 1 (/) + Kttf + 2/ b ) 

- «/ - 2/ 6 ) + 2K 2 (f) - K 2 (f + 2/ b ) + 2K 3 (f - 2/ 6 ) 

- 2JC 8 (/+2/ 6 )}, 
where 



*<» -£(-»" 53^ ©"*>/ 

w^^- 1 '" 275^(0""^' 
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The coefficients in the above equations are given as 
dn = Ci 2 , di 3 = 2cic 2 , du = c 2 2 + 2cic 3 
du = 2c 2 c 3 , die = c 3 2 
f?2i = c^, doo = 2c 4 c 5 , d 23 = c 5 2 + 2c 4 c 6 
r/04 = 2c 4 c 7 + 2c 5 c 6 , d 25 = c 6 2 + 2c 5 c 7 
d 26 = 2c 6 c 7 , dn = c 7 2 

$U = C1C4 , C?32 = C1C5 + C 2 C 4 , rf 3 3 = CiCe + C 2 C5 + C3C4 
«^34 = C1C7 + C 2 C 6 + C3C5 , d 3 6 = C 2 C 7 + C 3 C 6 , C?36 = C3C7 

013 = C1 3 , en = 3ci 2 c 2 , ei5 = 3cic 2 2 + 3ci 2 c 3 

0ie = 6cic 2 c 3 + c 2 3 , en = 3cic 3 2 + 3c 2 2 c 3 , ei 8 = 3c 2 c 3 2 

019 = c 3 3 , C21 = c 4 3 , 22 = 3c 4 2 c 6 

e 23 = 3c 4 c 5 2 + 3c 4 2 c 6 , 24 = 3c 4 2 c 7 + 6c 4 c 5 c 6 + c 6 3 

026 = 3c 4 c 8 2 + 6C4C6C7 + 3c6 2 C6 , e 2 6 = 6c 4 CeC7 + 3ceCe a + 3c B 2 C7 

027 = 3C 4 C7 2 + 6CbC 6 C7 + C 6 3 , 028 = 3c 5 C7 2 + 3C 6 2 C7 
029 = 3C 6 C 7 2 , 0210 = C7 3 , 032 = Cl 2 C 4 

C33 = 2cic 2 c 4 + Ci 2 c B , 034 = 02 2 04 + 2C1C3C4 + 2cic 2 c 5 + ci 2 c 6 
P35 = 2c 2 c 3 c 4 + C2 2 c 6 + 2cic 3 c 5 + 2cic 2 c 6 + Ci 2 C7 

036 = C3 2 C4 + 2c 2 c 3 cb + 2 2 C 6 + 2CiC 3 6 + 2CiC 2 C 7 

037 = Cz 2 C b + 2C 2 C 3 C 6 + Cxhn + 2CiC 3 C 7 

038 = 03 2 C6 + 2C 2 C 3 C7 , 039 = C3 2 C7 , 042 — CiC 4 2 

£ 43 = 2cic 4 c 5 + c 2 Ci 2 , C44 = C1C5 2 + 2cic 4 c 6 + 2c 2 c 4 c B + c z d 2 

045 = 2cic 4 c 7 + 2cic 5 c 6 + C2CB 2 + 2C3C4C5 + 2c 2 c 4 c 6 

046 = CiCe 2 + 2ciC 5 C7 + 2c 2 c 4 c 7 + 2c 2 c B c 6 + c 3 c 6 2 + 2c 3 C40e 
e 47 = 2ci0 6 07 + C2C6 2 + 2c 2 CbC 7 + 2c 3 04C 7 + 2c 3 c 5 c 6 
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^48 = C1C7 2 + 2ftftC 7 + ftft 2 + 2ftftC 7 

ftg = C2C7 2 + 2c 3 c 6 c 7 , e 4 io = ftc 7 2 
in = ft 2 , jiz = 2c 4 (c 5 - 2ci) 
•13 = (c 6 - 2ci) 2 + 2c 4 (c 6 - 4c 2 ) 
j 14 = 2c 4 (c 7 - 6c 3 ) + 2(c 5 - 2ci) (c« - 4c 2 ) 
Jib = (ft - 4c,) 2 + 2(c 5 - 2c x ) (c 7 - Gc 3 ) 
Jie = 2(c 6 - 4c 2 ) (c 7 - Gc 3 ), j i7 = (c 7 - 6ft) 2 
in = (ci + c 4 ) 2 , j'23 = 2(ci + c 4 ) (c 2 + 3c 6 ) 
Jm = (ft + 3c 6 ) 2 + 2(ci + c 4 ) (ft + 5ft) 
i 25 = 14c 7 (ci + C4) + 2(c 2 + 3ft) (ft + 5c 6 ) 
J26 = (ft + 5c 6 ) 2 + 14c 7 (c 2 + 3c 6 ), jn = 14c 7 (c 3 + 5c 6 ) 
ja = 49c 7 2 , j» = c 4 (ci + ft) 
Ju = ft(ft + 3c 5 ) + (ci + c 4 ) (c 6 — 2ci) 

j'33 = ft(ft + 5c 6 ) + (c 2 + 3c 6 ) (c 6 - 2ci) + (ci + c 4 ) (c 8 - 4c 2 ) 
; 34 = 7c 4 c 7 + (c 6 - 2ci) (c 3 + 5c 6 ) + (ft - 4c 2 ) (c 2 4- 3c 6 ) 

+ (c 7 - 6c 3 ) (ci + ft) 
J3B = 7c 7 (c 6 - 2ci) + (c 3 + 5ft) (c 6 - 4c 2 ) + (c 2 + 3c 5 ) (c 7 - 6ft) 
jM = 7ft (c 6 - 4c 2 ) + (c 3 + 5ft) (c 7 - 6c 3 ), j'37 = 7c 7 (c 7 - 6c 3 ) 
fcu = (ft + 2c 4 ) 2 , fci 3 = -2( Cl + 2c 4 ) (6ft - c 2 - 6c 5 ) 
Am = (6ft - ft - Oft) 2 - 2(ci + 2c 4 ) (20c 2 - ft - 10c 6 ) 
fcu = 2(6d - ft - 6ft) (20c 2 - c 3 - 10c 6 ) - 2( Cl + 2c 4 ) (42c 3 - 14c 7 ) 
ku = (20c 2 - c 3 - 10c 6 ) 2 + 2(6ci - c 2 - 6c B ) (42c 3 - 14c 7 ) 
ha = 2(20c 2 - ft - 10ft) (42c 3 - 14c 7 ), k m = (42c 3 - 14c 7 ) 2 
hi = c^, fc 22 = — 2c 4 (4ft — ft + 2c 4 ) 
fc 23 = (4ft - ft + 2c 4 ) 2 - 2c 4 (8c 2 + 12c B - ft) 
hi = 2(4c x - ft + 2c 4 ) (8c 2 + 12c 5 - ft) - 2c 4 (12c 3 + 30c 6 - c 7 ) 
/c 25 = (8c 2 + 12c 5 - ft) 2 - 112c 4 c 7 + 2(4ci - ft + 2c 4 ) (12c 3 + 30c 6 - c 7 ) 
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h 6 = 112c 7 (4ci - c 5 + 2c 4 ) + 2(8c 2 + 12c 5 - c 6 ) (12c 3 + 30c 6 - c 7 ) 

/c 27 = (12c 8 + 30c 6 - c 7 ) 2 + 112c 7 (8c 2 + 12c 5 - c 6 ) 

fc 28 = 112c 7 (12c 3 + 30c 6 - c 7 ), fc 29 = (56c 7 ) 2 

hi = c 4 (ci 4- 2c 4 ) 

fc 3 2 = — (ci + 2c 4 ) (4ci - c 5 + 2c 4 ) - c 4 (6ci — c 2 — 6c 5 ) 

/c 33 = -(ci + 2c 4 ) (8c 2 + 12c 5 - c 6 ) + (6ci - c 2 - 6c 5 ) (4ci - c 5 + 2c 4 ) 

- c 4 (20c 2 - c 3 - 10c«) 
A-34 = -(ci + 2c 4 ) (12c 3 + 30c 6 - c 7 ) + (6ci - c 2 - 6c B ) (8c 2 + 12c 5 - c 6 ) 

+ (20c 2 - c 3 - 10ce) (4ci - c 5 + 2c 4 ) - c 4 (42c 3 - 14c 7 ) 
k 35 = -56c 7 (ci + 2c 4 ) + (6ci - c 2 - 6c 5 ) (12c 3 4- 30c 6 - c 7 ) 
+ (20c 2 - c 3 - 10c 6 ) (8c 2 4- 12c 6 - c 6 ) 
+ (42c 3 - 14c 7 ) (4ci - c 5 + 2c 4 ) 
fc 36 = 56c 7 (6ci - c 2 - 6c 5 ) + (20c 2 - c 3 - 10c 6 ) (12c 8 + 30c 6 - c 7 ) 

4- (42c 3 - 14c 7 ) (8c 2 4- 12c 6 - c 6 ) 
A: 37 = 56c 7 (20c 2 - c 3 - 10c e ) + (42c 3 - 14c 7 ) (12c 3 + 30c 6 - c 7 ) 
fc 38 = 56c 7 (42c 3 - 14c 7 ), 
where 

ci = 2(a 2 / ( 2 + 2a 4 / b 4 + 3a 6 / 6 6 ) 
c 2 = -24(a 4 / 6 4 + 5a 6 f b 6 ), c 3 = 720a 6 / 6 6 
c 4 = a + az/ft 2 + flk/b 4 + ae/fc 6 
c 6 = -2(a 2 / 6 2 + 6a 4 /b 4 + 15a 6 / 6 6 ) 
c 6 = 24(a 4 / 6 4 + 15a 6 / 6 6 ), c 7 = -720a 6 / 6 6 . 
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